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>• ■ Recent observations 1 6 have revealed an unexpectedly high binary 

00 ' 

\ fraction among the Trans-Neptunian Objects (TNOs) that populate 

the Kuiper-belt. The discovered binaries have four characteristics they 
comprise a few percent of the TNOs, the mass ratio of their compo- 
nents is close to unity, their internal orbits are highly eccentric, and 
the orbits are more than 100 times wider than the primary's radius. In 
contrast, theories of binary asteroid formation tend to produce close, 
circular binaries. Therefore, a new approach is required to explain the 



unique characteristics of the TNO binaries. Two models have been 



proposed 7 8 . Both, however, require extreme assumptions on the size 
distribution of TNOs. Here we show a mechanism which is guaran- 
teed to produces binaries of the required type during the early TNO 
growth phase, based on only one plausible assumption, namely that 
initially TNOs were formed through gravitational instabilities 9 of the 
protoplanetary dust layer. 

The TNO binary with best known orbital elements, 1998WW31, has a mass ratio mijmx ~ 
0.7, eccentricity e ~ 0.8, and semi-major axis a ~ 2 x 10 4 km 2 . The inferred radii of 
the components are T\ ~ Wt% ~ 10 2 km, hence a/r% > 10 2 , in stark contrast to main 
belt asteroid binaries where m 2 /m! < 1, e ~ 0, and ajr\ < 10. Asteroid binaries in the 
main belt are probably formed by collisions 10 , as in the leading theory for the formation 
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of the Moon 11 - 12 , which naturally explains the extreme mass ratios and tight near-circular 
orbits 13 - 14 . 

In order to find a completely different process for TNO binaries, we can borrow ideas from 
dynamical binary formation in star clusters, where many roughly equal-mass binaries are 
found. There, we know of three formation paths: 1) two-body dissipational encounters 15 ; 2) 
three-body binary formation 16 ; and 3) exchange reactions 16 . Figure 1 schematically shows 
these paths. 

In the case of TNOs and asteroids, path 1 corresponds to the standard theory for binary 
formation, namely tidal disruption and giant impact during a close encounter (a and b in 
Figure 1). They indeed occur: each TNO has grown through accretion, and much of accretion 
has happened through collisions with an object comparable in mass to that of the growing 
TNO itself 17 - 18 . 

Path 2 (c in Figure 1) would require a near-simultaneous encounter of three massive objects 
with low enough velocities to allow an appreciable chance to leave two of the objects bound. 
For this to work, the random velocities of the most massive objects should be significantly 
lower than their Hill velocities, i>h = Rn^l(R), where Rn is the Hill radius of the asteroid 
at a distance R from the Sun and Q(R) is its orbital angular velocity. This path could 
work if there are ~ 10 5 100 km-sized objects embedded in a sea of small (< 1km) objects 8 . 
This assumption, however, is at odds with Goldreich and Ward's theory for the formation 
of planetesimals 9 through gravitational instability, and it is hard to see how objects in the 
Kuiper-belt could form from non-gravitational coagulation, because the time scales are far 
too long 19 . In contrast, the gravitational instability theory predicts the size of the initial 
bodies to be 10 — 100km. Starting with these bodies would make path 2 ineffective, because 
the velocity dispersion would be higher than the Hill velocity 17 - 20 . 

A combination of paths 1 and 2 has been proposed 7 . It was studied how a third massive 
body could capture a collision product of two massive bodies if the third body were near 
enough during the time of the collision. This mechanism seems unlikely to work, since it 
requires a number density of massive objects about two orders of magnitude higher than the 
value consistent with observations 8 . 

Another mechanism based on the dynamical friction from smaller bodies that can let a hyper- 
bolic encounter between two massive bodies produce a bound orbit has been proposed 8 . As 
we mentioned above, the gravitational instability theory for the formation of planetesimals 9 
would exclude the existence of such a sea of small objects. 

Path 3 (d in Figure 1) can operate on the binaries formed through path 1 (b in Figure 1), 
so we should check whether paths 1 and 3 together produce the right binaries in the right 
numbers. 

Consider a relatively massive TNO primary with mass m 1? in a binary with a secondary 
with mass m 2 <C m\. If the binary encounters another body with mass m 3 ~ mi, the 
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most likely result is an exchange reaction, in which the incoming object replaces the original 
secondary 21 . Figure 2 shows an example of such a reaction. 

If the incoming velocity of the third body is small, the binding energy of the resulting binary 
is comparable to or larger than that of the initial binary, since otherwise the secondary could 
not escape. In practice, the change in binding energy is relatively small, as will be shown 
later. Hence m^/flo ~ m\/a where a is the new semi-major axis after the exchange. This 
implies a/a « mi/m 2 3> 1. 

For the exchange reaction to occur, the third body should approach the binary to a distance 
comparable to a . The angular momentum carried away by the light body is small. Therefore, 
the periastron distance of the resulted binary is comparable to that of the initial parabolic 
orbit, in other words, to a . Thus, the eccentricity is e ~ 1 — a /a ~ 1 — m 2 /m 1 . 
We have run a series of numerical scattering experiments to obtain cross sections for different 
outcomes of binary-single body interactions (see Methods for details). The results show that 
80% of the cross sections corresponding to the change of the composition of the binary is 
accounted for by processes leading to the formation of binaries with two massive components. 
In figure 3 the distribution for the semi-major axis is strongly peaked around a ~ 20. This 
result is in good agreement with the simple argument presented above. This peak arises 
from the exchanges which take place at the first binary-single body encounter. The broad 
peak around a ~ 10 arises from "resonant encounters", in which encounters of two objects 
occur more than two times 23 . The eccentricity peaks at 0.95, as expected. In the case of 
1998WW31 with r\ = 75km, our length unit corresponds to 1500 km. In figures 2 and 3 we 
attach this physical scale to give a concrete idea of the size. Figure 4 clearly shows that the 
orbital elements of 1998WW31 are consistent with the binary having formed through the 
processes modeled here. 

We now confront our second task: to check whether the exchange path is efficient enough to 
produce the observed binaries. Starting with TNOs formed through gravitational instability 
of the protoplanetary dust layer 9 , the heaviest TNOs will accrete mass primarily through 
collisions with TNOs of comparable mass 17 - 18 . Many of these collisions are of the "giant 
impact" type, resulting in tight strongly unequal-mass binaries (b in Figure 1). Let us 
estimate how many such binaries are formed and what fraction of them are converted to 
wide binaries through exchange. 

We assume that one in three collisions between comparable TNOs gives rise to a binary 14 . 
The numerical simulations of planetesimal collisions yield a sizable formation rate of satellites 
of 43.3 % 14 . Since the relative velocity between objects in the TNO region is significantly 
lower than that assumed in their simulations, the corresponding binary formation rate in 
the TNO region is even higher. To stay on the conservative side, however, we have assumed 
a ratio of only 1/3. 

In 2/3 of the cases when no binary is produced, we have to wait for a typical time T 
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until another collision occurs. When a binary is formed, the cross section for subsequent 
interactions with a third body, taking into account gravitational focusing, is in proportion 
to a . Therefore, our newly-formed binary will undergo an exchange reaction on a time scale 
(r/a )T <C T, leading to a significant increase in a. Strong three-body interactions will 
subsequently occur on an shorter time scale (r/a)T <T. As a result, the semi-major axis 
will shrink systematically, with a 'thermal' distribution /(e) = 2e favoring high eccentricity 16 . 
When the orbit becomes small enough (r/a ~ 0.03), the change for collisions in resonant 
encounters becomes significant 24 . 

Let us assume that an exchange reaction turns a "giant impact" binary into a binary with 
a ~ 300r. Each subsequent strong encounter will on average shrink a by a factor 25 ~ 1.2. 
After a dozen encounters, a ~ 30r and a collision is likely to occur. The time scale for each 
encounter to occur is ~ (r/a)T, and the waiting time for the last encounter in this series is 
(1/30)T. Summing this series, we get a total waiting time of (T/30)/(l - (1/1.2)) = 0.2T 
before a collision. If three bodies collide, or if the third body escapes, the resulting system 
may be a single body (with an assumed chance of 2/3) or a strongly unequal-mass binary 
(a chance of 1/3). If the merger product and the third body are still bound, we have an 
equal-mass wide binary. 

Under these assumptions, in 1/3 of the cases, we wind up with an equal-mass TNO binary 
with the observed properties for a period ~ 0.2T, compared to a 2/3 chance to wind up with 
a single TNO for a period ~ T. Denoting by N s and N B the number of single bodies and 
the number of binaries, respectively, we have 



if we measure time in units of T. For the stationary state we have dN B /dt = dN s /dt = 0, 
and N B = 0.1N S . Therefore, the predicted binary fraction is ~ 10%. When accretion in the 
Kuiper-belt region diminished, the number of single and binary objects was frozen, with a 
ratio similar to this steady-state value. 

While our arguments are approximate, it is clear that after cessation of the accretion stage 
at least several percent of the TNOs were accidentally left in such a binary phase. The 
fact that more than 1% of the known TNOs are found to be in wide roughly equal- mass 
binaries is thus a natural consequence of any accretion model independent of the assumed 
parameters for the density and velocity dispersion of the protoplanetary disk or the duration 
of the accretion phase. We predict that future discoveries of TNO binaries will similarly 
show roughly equal masses, large separations, and high eccentricities. 
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Methods 



Numerical scattering experiments We performed 775,541 runs of numerical integrations 
of binary-single body encounters. We chose the relative velocity v ^ between the intruder 
(single body) and the binary at infinity to be much smaller than the orbital velocity of the 
binary components. In this regime, the scattering cross section a is expressed as a oc ctq/v^, 
because of gravitational focusing 7 - 25 . The parabolic orbit is the limiting case of — > 
where the nominal cross section a diverges, but o"o converges to a finite value. 
We fixed the initial orbital parameters of the binary, and changed the relative orbit of 
the intruder to obtain scattering cross sections. For encounters with finite v^, one should 
randomly change the direction of Voo and also randomly select impact parameter vectors uni- 
formly within a circle of given radius. This radius should be large enough so that encounters 
with impact parameter larger than this radius would not result in exchange or collision. For 
parabolic encounters, this procedure choice corresponds to choosing the periastron distance 
r p from a uniform distribution with an upper limit and selecting random orientations for the 
orbital plane and direction of the periastron. As an upper limit, we chose 20ao, where ao is 
the semi-major axis of the initial binary. Initial separations between the intruder and the 
binary were 100a . 

We use units in which G = m\ = m 3 = a = 1, where G is the gravitational constant, mi 
and 7723 are the masses of the heavier body of the initial binary and the intruder. The mass 
of the secondary component is rri2 = 0.05. The radii of the bodies are T\ = r 3 = 0.05 and 
i"2 = fi^/mi) 1 / 3 ~ 0.01842. The initial binary is circular. These values are typical for 
main-belt binary asteroids, with m 2 /mi < 0.1, and separations are 5 — 40 times the radius 
of the primary. 

We terminated the numerical integration when (a) one of the three bodies escaped beyond 
the Hill radius, or (b) when two bodies physically collided. For (a) we chose a Hill radius 
of 7"h = 100ao, since a typical value for a TNO with radius r, m — 10 21 g and R = 40 
AU is a = lOOOr, and the typical semi-major axis of a close, circular planetesimal binary is 
ao = 10r. When a collision occurred, the collision product was assumed to inherit the mass 
and motion of the center of mass of the two bodies. We adapted the radius so as to conserve 
the density. In both cases, we calculated the orbit of the remaining two bodies to determine 
the final state. 

We integrated the equations of motion of the three bodies in Cartesian coordinates. Inte- 
gration was done with an 8th order time-symmetric variable-stepsize Hermite method. To 
cross-check the results, we also performed two other sets of scattering experiments, one using 
the scattering package of START AB 22 , and the other using a 4th order Hermite scheme. In 
all cases, we used completely different programs for generating initial conditions, as well as 
for orbit integration and for the reduction of the results. All three experiments gave the 
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same outcome, within the expected fluctuations associated with the finite number of runs. 
There are six possible configuration changes that can take place as a result of a scattering 
event: (a) an exchange reaction resulting in a massive-massive binary; (b) an exchange 
reaction resulting in a massive-light binary; (c) a merger of bodies 1 and 2 resulting in 
a massive-massive binary; (d) a merger resulting in a twice-as-massive-light binary; (e) a 
merger of bodies 2 and 3 resulting in a massive-massive binary; (f) no binary is left, after 
three-body merging or two-body merging followed by escape. 

The cross section a for each event is as follows: (a) 12.1 (b) 1.3 (c) 0.9 (d) 1.3 (e) 0.9 (f) 
1.2 in our units. For our scattering experiments with r p < 20a , we found a cross section 
for preservation of the original binary of 1243.9 in our units, which implies that in 99 % of 
our experiments no configuration change took place. We also found that for r p > 20a no 
merger or exchange occurred, so our cross section measurements are complete. 
Another logical possibility would be: (g) all three objects become unbound. However, for 
that to occur the total energy would have to be positive. In our case, the low relative 
velocity of the incoming third body guarantees the total energy to be negative, hence (g) 
cannot occur. 

The physical reason for the dominance of (a) is the statistical tendency for equipartition of 
kinetic energy, which implies much higher velocities for much lighter bodies and therefore a 
much larger chance for those lighter bodies to escape. 

We acknowledge helpful comments on our manuscript by Peter Goldreich, Roman Rafikov, 
Re'em Sari, Keith S. Noll and Dan Durda. 
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Figure 1: Paths of formation of binaries. Path (a) is the formation through tidal disruption 
of one object followed by coagulation of fragments during a close encounter with the other; 

(b) a giant impact, where collision debris coagulates into a "moon" ; (c) an exchange reaction, 
where the incoming body replaces one component of a binary; (d) a combination of (b) and 

(c) . 
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Figure 2: An example of a binary-single-body exchange interaction. Bodies 1 and 2 having 
masses ni\ = 1 and ni2 = 0.1, respectively, form a binary with an initially circular orbit. 
Body 3, with mass = 1, encounters the binary on an initially parabolic orbit. In panel 
a, the whole scattering process is shown. Panel b shows the complex central interaction 
in more detail, while panel c shows the orbits of the initial (green) and final (red) binary, 
respectively. The final binary orbit is highly eccentric and much wider than the initial 
circular binary orbit. 
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Figure 3: Normalized differential cross sections for the formation of a 'massive-massive' 
binary. Normalized differential cross sections for the event-types (a), (c) and (e) in Methods 
are plotted with respect to the semi-major axis a (top panel), and eccentricity e (bottom 
panel) of the final binary. The physical units are given for reference at the top of the figure. 
The filled points are the total values, while the open circles are the contributions from the 
merger events (types (c) and (e) in Methods). 
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Figure 4: Distribution of orbital properties of 'massive-massive' binaries formed in our scat- 
tering experiments. Here a and e have the same meanings as in fig. 3. Contributions from 
exchange reactions, event-type (a) in Methods, are shown by green dots. They are limited 
by energy conservation to a < 20, and give rise to the horizontal rim in the middle of the 
figure. Contributions involving mergers, types (c, red dots) and (e, black dots), can lead to 
a values all the way to the Hill radius a fh 10 3 . The star symbol shows the observed orbit 
for 1998WW31. Boxes around the star indicate the observational 1- and 2-cr error bars. 
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